Eigenvalues and eigenfunctions of spin-weighted spheroidal harmonics 

in four and higher dimensions 



Emanuele BertQ 

McDonnell Center for the Space Sciences, Department of Physics, 
Washington University, St. Louis, Missouri 63130, USA 

Vitor CardoscQ 

Dept. of Physics and Astronomy, The University of Mississippi, University, MS 38677-1848, USA 

Marc Casalfi 

School of Mathematical Sciences, University College Dublin, Belfield, Dublin 4, Ireland 

(Dated: February 7, 2008) 

Spin-weighted spheroidal harmonics are useful in a variety of physical situations, including light 
scattering, nuclear modeling, signal processing, electromagnetic wave propagation, black hole pertur- 
bation theory in four and higher dimensions, quantum field theory in curved space-time and studies 
of D-branes. We first review analytic and numerical calculations of their eigenvalues and eigenfunc- 
tions in four dimensions, filling gaps in the existing literature when necessary. Then we compute 
the angular dependence of the spin-weighted spheroidal harmonics corresponding to slowly-damped 
quasinormal mode frequencies of the Kerr black hole, providing numerical tables and approximate 
formulas for their scalar products. Finally we present an exhaustive analytic and numerical study 
of scalar spheroidal harmonics in (n + 4) dimensions. 
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I. INTRODUCTION 



Spheroidal harmonics are special functions that play an important role in mathematical physics. The simplest 
prototype are the scalar spheroidal harmonics aSi m , obeying the differential equation 



[(1 - x 2 ) Si m , x ] x + 



(ex) 2 + Q Ai r , 



1-x 1 



oSlm — . 



(1.1) 



where l(= 0, 1, 2...) and m are integers, \m\ < I and oAi m is an angular eigenvalue to be found by imposing 
appropriate boundary conditions. The parameter c is, in general, complex. When c € R the eigenfunctions are 
called oblate; when c takes on pure imaginary values (c = ic/ with cj € IR) the eigenfunctions are called prolate 
Scalar spheroidal harmonics show up in a variety of physical situations ranging from light scattering |2( to nuclear 
modeling Jj|, signal processing and electromagnetic wave propagation There is a solid body of work on scalar 
spheroidal harmonics, some classic references being Stratton et al. Q , Meixner and Schafke and Flammer Q (see 
also 0, H, Ell for more recent developments on the subject). 

Generalized (four-dimensional) spin-weighted spheroidal harmonics (SWSHs) were first defined by Teukolsky [TlT | 
in the context of black hole physics. They result from the separation of angular variables in the equations describ- 
ing the propagation of a spin-s field in a rotating (Kerr) black hole background. Using the Kinnersley tetrad and 
Boyer-Lindquist coordinates, SWSHs satisfy a generalized version of Eq. (|l.ll) [see Eq. (|2.1|l below]. Now the an- 
gular separation constant, denoted by s Ai m , depends also on a spin-weight parameter s — 0, ±1/2, ±1, ±2 when we 
consider scalar, neutrino, electromagnetic or gravitational perturbations, respectively. When s — SWSHs reduce 
to ordinary, scalar spheroidal harmonics. When c = they reduce to spin- weighted spherical harmonics |12(, which 
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have innumerable applications in physics (see eg. |13| for an application to the analysis of anisotropies in the cosmic 
microwave background). The ordinary spherical harmonics are spin- weighted spherical harmonics with s = 0. 

Describing the angular dependence of a spin-s perturbing field in a Kerr black hole background, SWSHs and their 
eigenvalues find application in many theoretical studies of black hole physics. 

An important application concerns quasinormal modes of Kerr black holes 14]. The damped oscillation frequencies 
of Kerr black holes, u>, are uniquely determined by the black hole's mass M and specific angular momentum a. 
The calculation of quasinormal frequencies reduces to the solution of a coupled system of differential equations. 
One equation belongs to the class of generalized spheroidal wave equations, and describes the radial dependence of 
the perturbations; the other is a SWSH with c — auj, describing the angular dependence 0, 0, 0, 0, 0, l20| . 
For this reason, a detailed knowledge of the numerical and analytic properties of SWSHs is necessary to compute 
black hole quasinormal modes. Recent conjectures suggest a relation between highly damped quasinormal modes 
and black hole area quantization pll 12*2^. Verifying these conjectures for Kerr black holes calls for a calculation of 
the eigenvalues s A; m corresponding to values of c with imaginary part much larger than the real part, |c/| ^ |cr|. 
More phenomenologically, an important astrophysical problem is the determination of black hole parameters from 
gravitational wave observations [23l l24j . Gravitational waves emitted by the oscillations of a rotating black hole can 
be described by a superposition of quasinormal modes. The Laser Interferometer Space Antenna (LISA) has the 
potential to detect these waves with large signal-to- noise ratio (2j|. Given a detection, we would like to extract as 
much information as possible about the source. Ideally, we would like to determine M , a, the source location and 
the black hole's spin orientation from the observed waveform. An investigation of this issue requires the calculation 
of "scalar products" between different quasinormal modes, and in particular between the SWSHs describing their 
angular dependence [24[ . 

SWSHs also find application in quantum field theory in curved space-time. This theory is plagued with ultraviolet 
divergences: in particular, the expectation value of the stress-energy tensor, which is the crucial quantity that gen- 
erates space-time curvature via the semiclassical Einstein's equation, suffers from ultraviolet divergences and must 
be renormalized. This is usually a very difficult endeavour, particularly in the case of non-spherically symmetric 
space-times, which have recently been studied within this framework in 26] . An understanding of the high-frequency 
behaviour (both of the radial and angular parts) of the matter fields is essential in order to renormalize the various 
quantities in the theory |27|. 

Many theoretical scenarios rely on our Universe being (4 + n)-dimensional, with the n extra dimensions compactified 
on some small scale [2^]. The study of higher-dimensional black hole geometries, such as the Myers-Perry solution 
pg L requir es the introduction of higher- dimensional spheroidal harmonics (HSHs), that we shall define in Sec. IIIII 
|30ll3lL 132. l33| . TeV-scale gravity scenarios allow for the possibility that rotating, higher-dimensional mini-black holes 
are produced in particle accelerators such as CERN's Large Hadron Collider. These higher-dimensional, rotating 
black holes should evaporate by emission of Hawking radiation [30l l34| , and both SWSHs and HSHs are useful for a 
quantitative study of the evaporation process on the (3 + l)-brane and in the (4 + n)-bulk, respectively. HSHs also 
show up in recent studies of smooth geometries in the D1-D5 system [3^. Separating the Klein-Gordon equation 
leads again to five-dimensional spheroidal harmonics (even though they were not identified as such in |35|). 

There are no in-depth investigations of HSHs to date. On the other hand, the properties of SWSHs were (at least 
partially) investigated as soon as they were introduced. Press and Teukolsky [36J provided a polynomial fit in c of 
the eigenvalues s Ai m , which is valid up to c ~ 3. A formal perturbation expansion in powers of c was carried out by 
Fackerell and Crossman j3]j (see also j33, where some typos were corrected). An expansion for large real values of 
c has long been known for scalar spheroidal harmonics 0, @J3] • For general spin-s fields, the first attempt to find 
a large-c expansion was made by Breuer, Ryan and Waller |3iA l4p| . Their analysis, which was partially flawed and 
incomplete, has recently been revisited and corrected by Casals and Ottewill [23 • An expansion of the eigenvalue for 
large, pure- imaginary c is known for scalar spheroidal harmonics 0, an d was partially studied in [l9l l39| for 

general spin-s fields. However, as we already mentioned some of the studies for general spin-s were flawed, and we 
still lack a unified picture of the situation. 

The main purpose of this paper is to provide a complete analysis of general (four-dimensional) SWSHs, and to 
extend this understanding to HSHs. The plan of the paper is as follows. 

Sec. |n] contains analytic and numerical results for four-dimensional SWSHs. In Sec. Ill Al we define SWSHs and 
present a series solution first obtained by Leaver ^{|. In Sec. Ill Bl wc explain how this series solution provides us with 
a simple numerical algorithm to compute eigenvalues and eigenfunctions, and discuss alternative methods we used to 
check our results. In Sec. Ill Cl we give the series expansion of the eigenvalues for small c. In Sec. Ill Dl and lll El we present 
analytic results for large values of the argument in the oblate and prolate cases, respectively, supplementing these 
analytic expansions by numerical calculations. Sec. Ill Fl contains a numerical calculation of the prolate eigenfunctions. 
In Sec. Ill Gl we compute numerically the SWSHs at the slowly damped quasinormal frequencies of a Kerr black hole, 
providing numerical tables and a simple analytic approximation of their scalar products. While Sees. Ill Al III Bl III CI 
and III Dl are mostly review material, the results in Sees. Ill El III Fl and III Gl are new. 
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Sec. lIIIl is devoted to analytic and numerical calculations of the eigenvalues and eigenfunctions for HSHs. In Sec. lIII Al 
we review a series representation due to Ida et al. |3lj . illustrating the corresponding computational algorithm for 
the eigenvalues and eigenfunctions. The results in the following Sections are entirely new. In Sec. IIII Bl we give a 
series expansion of the eigenvalues for small c. Sees. IIII Ul and ITU Dl present series expansions for large values of the 
argument in the oblate and prolate cases, respectively, and back them up by numerical calculations. We conclude 
with a summary of relevant results and open problems. 

Other original results are contained in the Appendices. In Appendix^ we give an exact analytic solution of the 
HSH equation for c = 0, and use it to determine the number of zeros of the eigenfunctions in the region of physical 
interest. Appendix iBl provides a simple analytic solution of the HSH equation for special values of the parameters. 



II. FOUR-DIMENSIONAL SPIN- WEIGHTED SPHEROIDAL HARMONICS 



A. Series solution 



Using the Kinnersley tetrad and Boyer-Lindquist coordinates, the angular equation definining SWSHs results from 
the separation of the equations describing propagation of a spin-s field in the Kerr background fTT| : 



[(1 - X 2 ) s Slm tX ] x + 



(cx) — 2csx + s + s Ai m 



1 — X 



sS lm = 0, (2.1) 



where x = cos 9 and 9 is the Boyer-Lindquist polar angle. 

The angular separation constant s Ai m and the SWSHs s Si m are, in general, complex. They take on real values 
only in the oblate case (c € K) or, alternatively, in the prolate case (c = icj pure-imaginary) with s = 0. In the limit 
c — > the angular separation constant can be determined analytically: 

s A lm = 1(1 + 1) - s(s + 1) . (2.2) 

Some useful symmetry properties hold (see eg. [T^j): 

(i) Given eigenvalues for (say) positive to, those for negative to are readily obtained by complex conjugation: 

s Ai m = s A*_ m ; (2.3) 

(ii) Given eigenvalues for (say) negative s, those for positive s are given by 

- s A lm = s A im + 2s. (2.4) 

Exploiting these symmetries, in our numerical calculations we only consider s < and m > 0. 

(hi) Let us define p = ic. If p and - s Ai m corresponds to a solution for given (s, I, to), then another solution can 
be obtained by the following replacements: to — > — m, p — > p* , _ s A; m — >_ s A*_ m . 

(iv) In the prolate case, if s Si m is a solution with eigenvalue s Ai m and cj > 0, then s Sf m is a solution with eigenvalue 
s A* lm and ci < 0. 

Leaver found the following series solution for the angular eigenfunctions [l5j : 

oo 

s Si, n (x) = e cx (1 + x) k - (1 - x) k+ Op(l + xf , (2.5) 

where k± = \m ± s\/2. The expansion coefficients a p are obtained from the three term recursion relation 

a Q ai + P ao = 0, (2.6) 
pap+i + f3 p a p + 7 P a p _i =0, p = l,2... (2.7) 



a 



with 



a p = -2(p + l)(p + 2A- + l), (2. 

/3 P = Pip - 1) + 2p(fc_ + fc+ + 1 - 2c) 

- [2c (2A:_ + s + 1) - (*_ + fc+) + k+ + 1)] - [c 2 + s(s + 1) + s A ;m ] , 

7 P = 2c (p + k- + k+ + s) . 
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Given a (generally complex) argument c, the separation constant s Ai m can be obtained solving numerically the 
continued fraction equation 



a 7i ai7 2 « 2 73 _ n , n . 



or any of its inversions . 



B. Numerical calculation of the eigenvalues and eigenfunctions 

Flammer's classical reference pj is largely dedicated to tabulating eigenvalues and eigenfunctions of scalar spheroidal 
harmonics. The main limitation of Flammer's impressive work is that his tables only deal with the prolate (pure 
imaginary frequency) and oblate (pure real frequency) cases. The general case of complex frequencies (which is of 
interest, for example, in the calculation of the angular dependence of quasinormal modes) is not covered. A pioneering 
numerical work in this sense was carried out by Oguchi 0, who computed angular eigenvalues for complex values 
of c and s = 0. Quite recently this topic received more attention. Li et al. published a useful review of numerical 
methods to compute eigenvalues and eigenfunctions for s = ,§]. Falloon et al. developed a Mathematica notebook 
to compute s = harmonics for general complex values of the frequency 9] . Finally, Barrowes et al. provided a rapid 
and accurate method to calculate the prolate and oblate scalar spheroidal wave functions and their eigenvalues for 
complex frequencies in the limit |c| — > oo |ldj| . 

In comparison, numerical calculations of the eigenfunctions and eigenvalues for general spin s are scarce. To make 
things worse, until recently the few analytic predictions were contradictory (the situation for large real frequencies 
has finally been clarified in |27| . while large pure imaginary frequencies were considered in [T^L 

Leaver's solution gives a simple and practical algorithm for the numerical calculation of eigenvalues and eigenfunc- 
tions. Start from the known analytic eigenvalue for c = 0, Eq. I|2.2[) . Use this as an initial guess, increase the value of c 
and solve numerically Eq. I|2.9[) to get the eigenvalue for c ^ 0. Once the eigenvalue is known, compute any number of 
series coefficients a p using the recursion relation, and plug them into the series solution (|2.5|) to get the corresponding 
eigenfunction to any required precision. In our numerical calculations we truncate the series at p = p m ax = 1000 (but 
usually an excellent approximation to the eigenfunction can be obtained keeping only ~ 10 terms). 

This algorithm only determines the eigenfunction up to a normalization constant, which can easily be fixed by 
imposing the normalization condition 

J \ s S lm (x)\ 2 dx = l. (2.10) 

We independently checked our numerical results for the eigenvalues and eigenfunctions using different methods. 
These methods are variants of those described in detail in [27j, where they were used to deal with oblate SWSHs, 
so we only describe them briefly. To check the eigenvalues from Leaver's continued fractions we adapt the procedure 
derived by Sasaki and Nakamura ^l| for s = — 2 and c S R. The idea is to approximate the angular differential 
equation by a difference equation, which can be written in matrix form. The eigenvalue is then obtained by imposing 
the determinant of this (tridiagonal) matrix to be zero. The Sasaki-Nakamura method faces a numerical difficulty as 
the discretization grid is refined: the value of the determinant becomes very large, making it increasingly difficult to 
obtain the eigenfunctions. Therefore, in order to find the eigenfunctions we complement the Sasaki-Nakamura analysis 
by an adaptation of the shooting method introduced in |42j for the special case of (scalar) spheroidal harmonics with 
real frequency. We start the integration of the angular equation with arbitrary values for (i) the eigenvalue, (ii) the 
angular function at x = —1. The eigenvalue results from imposing the derivative of s Si m near x = +1, as found by 
numerical integration, to agree with a certain analytic approximation of the derivative, coming from a power series 
expansion around x = +1. In |27|. for large real frequency it was sometimes necessary to look for an extreme (rather 
than a zero) of either the determinant or the difference of derivatives near x = +1. That is because in the oblate 
case the eigenvalues for some modes "pair up" for large frequency. This degeneracy does not occur for large and 
pure-imaginary frequency (prolate case). On the downside, when using the "shooting" method we must now find the 
zero of a complex function. We locate this zero using Newton's method fyj, which is globally convergent. 

As a first example of the application of Leaver's method, in Fig. ^ we show some prolate and oblate eigenvalues for 
s = 0. Progressively increasing the inversion index of the continued fraction |l5j | we can track the eigenvalues (which 
in this special case are real) up to quite large values of c. Our numerical code passed a number of tests. We first 
verified that our results are in agreement with Tables 10-12 of in the prolate case and Tables 130-132 of in the 
oblate case. For complex c we were able to reproduce Table 2 in JuL As a final sanity check, we verified in a few 
representative cases that the Mathematica notebooks presented in |8|, |9j are also in agreement with our code. 
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FIG. 1: Angular eigenvalue oAi m as a function of a (prolate case) or c (oblate case), for selected values of I and m. 



Notice that prolate and oblate eigenvalues coincide when c — * 0, but their behavior for large frequencies is markedly 
different. Prolate eigenvalues grow linearly with |c/| for large |cj|, while oblate eigenvalues are proportional to — c 2 
for large c. In the following Section we derive analytically asymptotic expansions in the limit of small c, large real c 
and large pure-imaginary c. 



C. Small-c expansion 



The analytic properties of the ei genv alues and eigenfunctions of the SWSH equation (|2.1I) have been studied by 
many authors Q, 0, |27], H(| 1371 l38l l39l Eof . Here and in the following Sections we present for the first time a 
compact and complete summary of the analytic results for the eigenvalues in various asymptotic regimes. 

The small-c expansion was worked out many years ago for the scalar case s = (se e eg. |l| ) . The general s case was 
considered by Press and Teukolsky [3(|, Fackerell and Crossman [37| and Seidel [3^|. Here we present an alternative 
and simpler derivation of their results, which can easily be extended to the higher-dimensional case (see Sec. IIIIB)l . 
Our method is directly based on the continued fraction representation. For c = 0, the eigenvalue s Ai m is explicitly 
determined from the requirement that the series expansion has a finite number of terms, since otherwise it is divergent 
[IH . Taking c = in \2.9\ and imposing (3 r — for some integer r > we obtain 

s A lm = (r + H)(r + |m| + 1) - s(s + 1) . (2.11) 

Appendix Ogives an alternative derivation of this result in a more general framework. Setting r = I — |m| we recover 
the result (|2.2|l : furthermore, the condition r > implies the well-known constraint on the angular quantum numbers, 
l>\m\. 

For finite but small c one can use standard perturbation theory to express the separation constant as a power series 
in c. The standard procedure is somewhat involved |37ll38|. but it can be greatly simplified as follows. For c = we 
know the recursion relation l|2.6|l has a finite number of terms r. Therefore, when expanding around c = we find it 
convenient (although strictly not necessary) to use the r-th inversion of Eq. I|2.9|l : 

P _ O-r-llr QV-27r-l Qp7l _ a rlr+l CX-r+lJr+2 ^ ^ 

Pr-1 Pr-2— Po Pr+1~ Pr+2 — 

Now we expand the separation constant as a Taylor series: 

oo 
p=0 
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where /q = 1(1 + 1) — s(s + 1). Plugging l|2.13|) into Eq. (|2.12|l . expanding the resulting expression in powers of c and 
defining 

1 i uis 

-(a + 0) =max( H ,N), -(a-/3) = ^ ffly , (2.14) 

we have for the first coefficients 

fo = + + (2.16a) 

* - -WTVy < 2 - 16b ' 

h = h(l+l)-h(l)-l, (2.16c) 

2fe(Z)ms 2 2h(l + l)ms 2 

/3 = 71 TT7271 I il ~ Tn i Twi I o\ ' (2.16d) 



(Z — 1)Z 2 (/ + 1) /(/ + l) 2 (Z + 2) 
, 4 / 4h(l + l) 

^ 2 G + l) 4 (Z + 2)2 (/- 

/i 2 (/ + l) h(l)h(l + l) h 2 {l) (I - l)h(l - l)h(l) 
2(1 + 1) + 2l 2 + 21 2~T + 4/ 2 - 21 



2 4 / 4h(l + l) 4h(l) \ (l + 2)h(l + l)h(l + 2) 

h ~ 171 3 \p(i + iy(i + 2)i (i-iyi^i + iy) 2(/ + i)(2z + 3) + 



36 f 8h(l) 8h(l + l) 

f 5 = ms 



ms h(l) 



l 6 (l + l) 3 (l-l) 3 P(l + l) 6 (l + 2)\ 

h(l + l)(7l 2 + 71 + 4) h(l-l)(3l-4) 



l 3 (l + 2)(l + l) 3 (l - 1) l 3 (l + l)(2l - l)(l - 2) 

2 



/ (3/ + 7)h(l + l)h(l + 2) 3h 2 (l + 1) 3h 2 (l) 



h = 



\l(l + l) 3 (l + i)(2l + i) 1(1 + If (l + 2) P(l-l)(l + l) 
16m 4 s 8 / h(l + l) h(l) 



l 4 (l + l) 4 \(l + If (I + 2) 4 l 4 (l - 1) 4 / 
4m 2 s 4 / (3Z 2 + 14Z + 17)/ i (Z + 1)^(Z + 2) ih 2 (l + l) 3h 2 (l) 



(2.16e) 



(2.16f) 



l 2 (l + l) 2 \ (l + lf(l + 2)(l + 3) 2 (2l + 3) (l + lf(l + 2f l 3 (l-lf 

Am 2 s 4 Hill 4 + 221 3 + 311 2 + 201 + 6)h(l)h(l + 1) (3l 2 - 81 + 6)h(l - l)h(l) \ 
l 2 (l + l) 2 \ l 3 (l-l) 2 (l + l) 3 (l + 2) 2 + Z 3 (Z — 2) 2 (Z — 1)(2/ — 1) ) + 

( i± ^ + £(" + w + " " <* + 10 >"« + " + ) + 

- ^ + mw + « + 3 »""» - ^ + « + » - ( ' 2 " ,(3 ';^ 2) """' ) + 

^fciM ( (i _ 1)(n _ 3)M0 _ (1 _ _ 1} _ , (2 , 6g) 

in agreement with |3Sj . Notice that our separation constant s Ai m differs from the separation constant s Ei m of 
Refs. |H|3|: s A lm = s E lm - s(s + 1). 



D. Large and real c (oblate case) 



For s = 0, analytic power-series expansions for large (pure real and pure imaginary) values of c have long been 
known P, [|| . In this Subsection and in the next we review and extend known analytic results for general values of 
the spin- weight s. Large-c expansions of the eigenvalues are usually based on the following idea. First, an asymptotic 
approximation of the solution is found. An asymptotic expansion of the eigenvalue is then obtained by matching the 
number of zeros of the asymptotic solution with the (known) number of zeros of the exact solution (as determined, 
eg., in pfljl'). This matching relies on the solution being analytic in c. 
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Breuer, Ryan and Waller jijj perform an asymptotic analysis for large, real c (and fixed to) of the SWSH equation. 
Their analysis is based on similar analyses in the s = limit Q, ||J and on a partially-flawed previous study for general 
spin |39j |. It yields the following asymptotic expansion for the separation constant: 



sAi m = -(? + 2 s q lm c - i [ s qf m - to 2 + 2s + ll + -Ai~ 
2 c 



1 



[5 s qf m - (6to 2 - 10) s qf m + to 4 - 2to 2 - 4s 2 ( 8 qf m - m 2 - l) + l] 



where 



64c 2 
\a 3 + O (1/c 4 



[ s qf m - m 2 s qi m + s q lm ~ s 2 ( s q hn + to)] 

( s qim - m - 1 + 2s) ( s qi m - to - 1 - 2s) ( s qi m - m - 2s - 3) x 
x { s qi m - to + 2s - 3) ( s qi m + to - 3) 2 ( s qi m + to - l) 2 - 

- { s qim - to - 2s + 1) ( s qi m - TO + 2s + 1) ( s q im - m - 2s + 3) { s qi m - m + 2s + 3) X 
x ( s qi m + to + l) 2 ( s ^ m + TO + 3) 

{sQlm - m + 2s - 1) ( s q ;m - m - 2s - 1) ( a q lm - l) 2 ( s gf im + to - 1)^ 
(sQlm + m + l) 2 { s qi m + if { s qim - m - 2s + 1) ( s q lm - to + 2s + 1) 



(2.17) 




(2.18) 



2A X 



{aQlm - to + 2s - 1) ( s qi m to 2s 1) ( s g; m + m - 1) + 



+ ( s tR m + to + l) 2 ( s qi m - to - 2s + 1) ( s gv - TO + 2s + 1) 

Unfortunately, Ref. 01 left the parameter s qi m undetermined. The recent work by Casals and Ottewill [27} 
determines s qi m through a complete asymptotic analysis for large, real c (and fixed to). First they find an asymptotic 
expansion for the angular solution valid everywhere in x G [— 1,+1], Then they determine s qi m by imposing the 
asymptotic angular solution to have the right number of zeros. This could not be done successfully by Breuer, Ryan 
and Waller as their asymptotic angular solution was not valid near the origin x = 0, where the solution may have one 
zero. The value of the parameter s qi m obtained in |27| is: 



sqir. 



I — max(|m|, |s|) + 



(|to 



+ l-z 



s qi m = 2(1- max(|m|, |s|)) + \m qp s\ T s + 1 
where s l m = max(|TO|, |s|) + (\m + s\ — \m — s|)/2 + s, and the variable 



if I > max(J m , - s l m ) 
if I < ± s lm 



if I -J, 

1 if I - J r 



even 
odd 



(2.19a) 
(2.19b) 

(2.20) 



gives the number of zeros of the asymptotic angular solution at x = 0. 

In Fig. |21 we use Leaver's method to give numerical support to the analytic results, computing selected oblate 
eigenvalues for s ^ 0. Fig. [21 only shows the lowest radiatable multipoles (that is, the lowest values of I such that 
I ^ l s l)i which are physically the most relevant. Numerical results fully confirm the analytic expectations. The 
imaginary part of the eigenvalues is always zero within the numerical accuracy. The small-c and large-c behaviors are 
both in agreement with the analytic predictions. As an additional code check we also verified that the numerically 
computed eigenvalues satisfy the symmetry relation H2.4(l . 



E. Large and pure-imaginary c (prolate case) 



In the scalar case s = 0, the asymptotic analysis of the eigenvalues for large and pure-imaginary c presents no 
particular difficulties, and can be found in A very accurate approximation of the numerical eigenvalues in the 
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FIG. 2: Angular eigenvalue for the lowest radiatable multipoles with s = — 1 (left) and s — —2 (right) in the oblate case. 



left panel of Fig.^ which is valid within a few percent down to c ~ 2, is: 



A lm = (2L + l)|c/| - (2L 2 + 2L + 3 - 4m 2 )2" 2 - — (2L + l)(L 2 + L + 3- 8to 2 )2~ 

\ci\ 

:[5(i 4 + 2L 3 + 7L + 3) - 48m 2 (2L 2 + 2L + 1)]2~ 6 + + C(1/| C/ | 3 ) , 



(2.21) 



where L = I — \m\ is the number of zeros of the scalar harmonics in the interval — 1 < x < 1. An expression accurate 
to 0{l/\c,\ 6 ) can be found in Eq. (8.1.11) of 0. 

Ref. (33] provides, to our knowledge, the first attempt to generalize Eq. (|2.21|) to general spin s. Unfortunately the 
number of zeros of the SWSH given on page 115 of]39j is wrong, and was later corrected by Breuer et al. (Theorem 
4.1 in 40]). Therefore, although the derivation in 1391 is essentially right, the expression for the leading term is not. 

Breuer's analysis for general spin was corrected in [13], but two open issues remain. In the scalar case, the asymptotic 
solution found in 0, [|| is a valid approximation to the spheroidal harmonics only in a region far from the end-points 
x = ±1. We expect the function, which is real, to have most zeros far from the end-points; still, in principle we 
cannot discard the possibility that the function does have a zero near the end-points. As mentioned in the beginning 
of Sec. Ill Dl the omission of a zero would lead to the wrong asymptotic expansion for the eigenvalue (this is indeed 
what happened in (Ijj m the case of large, real c and general spin s). 

For general spin not o nly we run into the same problem, but we have the additional complication that the solution 
is complex. However, in |!9ll39j] the number of zeros of the exact, real solution for real c was matched with the number 
of zeros of the real part of the asymptotic (complex) solution for pure-imaginary c. As indicated in P, Q for the case 
s = 0, the number of zeros of the SWSH is the same for c 2 € (— oo,+cx}), but this is only true for real functions. 
Numerical calculations (see Sec. Ill Fl below. in particular Figs. |5] and |HJ show that the number of zeros of the real part 
and of the modulus of the SWSHs are both functions of c (when c is pure-imaginary) . Therefore it is not completely 
clear why the asymptotic analysis in 19, 39] for the eigenvalue should be correct. 

Below we provide analytic and numerical arguments in support of the result obtained in |l9l I39I ] . We derive an 
asymptotic solution for general spin-weight s which is a valid approximation to the SWSH near the end-points. This 
"outer" solution does not contain any zeros in the case s — 0, and this explains why the asymptotic behaviour of the 
eigenvalue is correct, at least for scalar perturbations. We could not generalize this argument to electromagnetic and 
gravitational perturbations, but we give numerical evidence that, in all cases we considered, the asymptotic behaviour 
of the eigenvalue given in [19] for \s\ = 1 and |s| = 2 is correct. 

In order to find an asymptotic solution we define a new angular wavefunction s yi m (x) via 



sSi m (x) = (1 - xf+ (1 + x) k - s yi m (x) , 



(2.22) 
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and change the independent variable setting u = \j2\ci\x. Substituting into Eq. (|2.1() we get 



(2\a\ -u 2 )^-2 [^2M (*+ - *-) + ( k + +*- + !) '■ 



d 

du 

(2.23) 

\ C r\ U 2 I 

+ s A lrn + s(s + 1) - (k+ + fc_) + + fc_ + 1) - J-i iv/2Msu \ s y lrn = 

When |c/| — > oo and in the region l/|c/| <C |a;| <C 1, the equation satisfied by the real part of s yi m becomes a 
parabolic cylinder equation. Therefore the real part of the "inner solution" s yf^, which is a valid approximation in 



this inner region, is a parabolic cylinder function Dr, with L zeros: 5ft ( s y™n) ( x ) = ®L ^^/2|c/|a;J , where 

s A lm = {2L + l)\ Cl \ + O(\a\ ) , |cj| -> oo . (2.24) 



This is essentially how the inner solution was found in [lj, |6j, [Hi |39j . Note, however, that the slightly different change 
of variable (|2.22|) means that the inner solution s S l [ ^ 1 (x) we have obtained is always regular at x = ±1, whereas the 
one given in 19, 39] is not. 

It is known 43] that the zeros of the parabolic cylinder function D L (*/2\ci\x^ lie within |x| < \J (2L + l)/|c/|, so 

they are all within the inner region. We can determine the value of L by equating the number of zeros of the parabolic 
cylinder function to the number of zeros of the real part of the exact solution for pure-imaginary c — ic/ (i.e., the real 
part of the prolate SWSH) in the region l/|c/| <C |x| <C 1. Unfortunately, the latter number is not known. In [l9ll39| 
the number of zeros of the parabolic cylinder function is instead equated to the number of zeros of the exact solution 
for real c (i.e., the oblate SWSH) in the whole region x E (—1, +1), as given in [4fJ. The value of L so determined 
turns out to be 

L = Z-max(|m|,|s|). (2.25) 

Higher order corrections in the asymptotic expansion can be obtained as indicated in [jj. We checked Eqs. H2.24J1 and 
(12.25(1 numerically, solving Leaver's angular continued fraction for s Ai m as a function of the complex parameter c. 

Now we take one more step towards understanding the asymptotic behaviour given by Eqs. (|2.24|l and (|2.25|) 
providing an "outer" asymptotic solution, which is a valid approximation to the SWSH in a region far from the origin. 
A WKB-type expansion similar to the one in [22j gives the result: 



s S%?(x) = s aim (1 - x 2 ) 1/4 x L (l + Vl- x 2 ) L V2 (x - iVl - z 2 ) V^I^V 
Mm (1 - XT' 7 ' x-^ 1 (l + VT^) L+1/2 (x iVl - ^) +S e-l-lv^ 



(2.26) 



where s ai m and s bim are constants of integration. This solution is valid when |i| » l/y/\cj\ and |x| ~ 1 — 0(|c/| e ), 
with —1 < e < 0. Matching the outer and inner solutions so that they coincide in their overlap region, we get the 
following final expression for the outer solution (valid near x = ±1): 

s5 ™ t ' ±1 (x) = H)*2 L+1 /2 (±V2M) L e- 3 ^ 2 (1 - x 2 Y 1,A x L (l + VT^y L ^ /2 (x-\^/T^y S e+\ c '\^^ . 

(2.27) 

This solution formally diverges at x = ±1, and actually it is not valid at the end-points, but only in a region 
arbitrarily close to the end points. In such a region the exponential decay with \ci\ takes care of the diverging factor 

(l — a; 2 ) 1 ^ 4 . It can easily be checked that: (1) the outer solution has no zeros when s — 0; (2) its real part has no 
zeros within x £ (—1, +1) when |s| = 1; (3) its real part has two zeros located at x = ±l/y/2 when |s| = 2. 

We finally have a complete account of scalar perturbations (s = 0). The outer solution has no zeros, so all zeros of 
the asymptotic solution are provided by the inner solution. Equating the number of zeros of the parabolic cylinder 
function to that of the oblate SWSH is therefore valid in this case, which is why Eqs. I|2.24|l and H2.25fl are correct 
for s — 0. For s ^ we are one step closer to understanding the asymptotic behaviour. Now we know all the zeros of 
the real part of the asymptotic solution, and we "only" need to find out the number of zeros of the real part of the 
SWSH for pure-imaginary frequency. We have numerically checked in many cases that Eqs. I|2.24|l and l|2.25|l give 
the right behaviour, and we believe them to be correct even for s ^ 0, but only an analytic proof can settle the issue 
once and for all. 
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FIG. 3: Prolate angular eigenvalues for s = — 1 (left) and s = —2 (right). At variance with Fig. [21 the angular eigenvalue a Ai m 
is now complex. Lines limiting to s Ai m = 1(1 + 1) — s(s + 1) as \c\ — > are the real parts of the eigenvalues. Lines approaching 
zero as |c| — » are the imaginary parts if ci < (or their modulus if c/ > 0: see property (iv) of Sec III A0 . 



Fig. |31 shows prolate eigenvalues s Ai m for s ^ and a few selected values of (l,m), providing numerical support 
to our analytic conclusions. The angular eigenvalue s Ai m is now complex (unless m = 0). We can still limit our 
calculations to positive m's because of the symmetry property l|2.3[l . From Fig.|31we see that in the prolate case the 
real part dominates as c/ — > oo. The numerical results for any s are consistent with the linear growth predicted by 
Eqs. (|2.24|) and H2.25fl . A peculiar feature can be seen in the real part of -2^22, which seems to "bend down" at 
C] ~ 3. We verified that the presence of this "bending" is not a numerical artifact using both Leaver's method and 
the shooting method. There is no discontinuity at the bending location, just a smooth change of slope. Ref. |Tsj 
found that, in the intermediate damping regime, the real part of the Kerr quasinormal frequencies with I — m = 2 
tends to lor = 2Q, Q being the angular velocity of the horizon, unlike the frequency of modes with any other value of 
I and m. We suspect this anomalous behavior of quasinormal modes with I = m = 2 and s = —2 in the intermediate 
damping regime could be related with the observed bending of the angular eigenvalue. 




FIG. 4: K_2^422 as a function of |c| for different, fixed values of the phase angle 6. 

Finally, Fig.0]is a computational exercise to illustrate in a specific case the transition from the prolate to the oblate 
regime. We set c = |c|e and compute the real part of -2^22 along trajectories in the complex plane with fixed 



11 



values of 9 G [0, tt/2]. A phase angle 9 — corresponds to the oblate case (Fig. 0) and a phase angle 9 = 7r/2 to the 
prolate case (Fig.[3J). For intermediate values of 9 the function 3?_2^422(|c|) smoothly deforms from the prolate to the 
oblate character. Notice also that the bending of the I — m = 2 prolate eigenvalue we see at \c\ ~ 3 (see also Fig.^J 
disappears immediately as 9 becomes nonzero. In our opinion, this provides some circumstantial evidence that the 
bending is not due to a branch cut. We refer to Q, E3 for an extensive discussion of branch cuts and their effect on 
the prolate/oblate nature of the eigenvalues when s = 0. 

F. Numerical calculation of the prolate eigenfunctions 

Once we have obtained the eigenvalues it is a simple matter to compute the eigenfunctions for given values of s, 
I, m and c: we just need to compute the coefficients a p from the recursion relation (|2.6|l - l|2.8|l and plug them into 
Leaver's series solution, Eq. (|2.5|) . The eigenfunctions for s — are well known pj, and so are the eigenfunctions 
for oblate harmonics (real c) with general spin s |27|. Therefore here we concentrate on prolate eigenfunctions with 
different values of cj. In black hole perturbation theory, according to Leaver's convention on the Fourier transform 
|15| stable perturbations correspond to quasinormal frequencies with ljj < 0. For this reason, below we pick cj < 0. 
Using the symmetry property (iv) of Sec. Ill Al eigenvalues and eigenfunctions for c/ > are trivially obtained by 
complex conjugation. 
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FIG. 5: Modulus of the eigenfunction for prolate SWSHs with s = 
the bottom row to I = 3. Left to right: cj = 0, — 1 and —5. 



-2 and different values of m. The top row refers to I = 2, 



The modulus of prolate eigenfunctions with increasing values of |c/| and s = —2 is plotted in Fig. [S] The top 
row refers to I = 2, the bottom row to I = 3. Each panel shows eigenfunctions for different (positive) values of m. 
The left panels refer to the simple case ci = 0. In this limit the eigenfunctions reduce to spin-weighted (s = —2) 
spherical harmonics, which can be thought of as being prolate and oblate at the same time [l^. The "spikes" in this 
logarithmic plot are zeros of the eigenfunctions, and their number is in agreement with Eq. I|2.25|l for cj = 0. However, 
as anticipated, Fig. 0shows that (for s ^ 0) as soon as cj ^ the number of zeros in the region —1 < x < 1 changes. 

Fig. El only gives the modulus of the eigenfunctions, to illustrate the point that prolate harmonics with spin- weight 
s = —2 have no zero for — 1 < x < 1 when cj 0. In Fig. [S]we show both the real and imaginary part of prolate 
eigenfunctions for cj = —5. Their modulus can be read off the top right panel of Fig. [S] The number of zeros of 
the real part of the SWSH is not zero, as predicted by Eq. I|2.25|l for \m\ < |s| = 2 and I = 2. This illustrates that 
the number of zeros of the real part of the eigenfunction (as well as the number of zeros of its modulus: see Fig. EJ 
depends on c/. 
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FIG. 6: Real (left) and imaginary part (right) of prolate I = 2, s = —2 eigenfunctions for a — —5 and m > 0. Note that the 
number of zeros of the real (or imaginary) part is not zero, as predicted by Eq. (I2.25H : the number of zeros of the real part 
varies with c/. 



G. Scalar products of the eigenfunctions at the Kerr quasinormal frequencies 

As a useful application to general relativity of the numerical algorithm described in Sec. Ill Bl here we study the 
angular eigenfunctions corresponding to different quasinormal mode frequencies and different values of a/M, for 1 = 2. 
To our knowledge, these eigenfunctions have never been shown before in the literature. The procedure we use is by 
now standard and can be found in many papers 0, 0, 0, 0, |2(| , so here we only sketch the basic idea. In the Kerr 
space-time, linear gravitational perturbations are described by a pair of coupled differential equations: one for the 
angular part of the perturbations, and the other for the radial part. The radial equation resulting from the separation 
of the field equations in the Kerr background is given (eg.) in |lllll5j. The angular equation is a SWSH equation with 
c = auj. Boundary conditions for the two equations can be cast as a couple of three-term continued fraction relations 
of the form (|2.9I) . Computing quasinormal frequencies u) and the related angular eigenvalues is simple. For assigned 
values of s,l,m,a and u>, first find the angular separation constant s A[ m looking for zeros of the angular continued 
fraction. Then replace the corresponding eigenvalue into the radial continued fraction, and look for its zeros as a 
function of to. 

Once we know the quasinormal mode frequencies and the angular eigenvalues, we can find the angular eigenfunction 
following the procedure of Sec. Ill Bl The results are plotted in Fig.|7|for different gravitational (s = —2) quasinormal 
frequencies with I = 2, |m| < I and different values of a/M. The modulus of the eigenfunction is scarcely affected 
by rotation and almost unaffected by the overtone index, at least for slowly damped modes. The phase (hence the 
relative magnitude of the imaginary part of the eigenfunctions) is more sensitive to angular momentum and mode 
damping. It is generally larger for counterrotating modes, and it grows with the overtone index N. In general, the 
magnitude of the phase seems to be related with the magnitude of the imaginary part of the quasinormal frequencies. 
For corotating modes in the extremal limit the imaginary part of quasinormal frequencies goes to zero, as first shown 
by Detweiler |44| . This explains what, at first sight, might look like a puzzling feature of the plots: for near-extremal 
Kerr black holes (a/M = 0.98, green lines) the phase is actually smaller than for black holes with a/M = 0.80 if the 
modes are corotating. The opposite is true for counterrotating modes, whose imaginary part does not tend to zero in 
the extremal limit. 

In analysing data from gravitational waveforms emitted by an oscillating black hole it is necessary to compute the 
scalar products [24| 

J -2<SYm* (a^lmN) -2'S'/'m' (a07/' m ' N')dfl — a m W NN> {a)&m,m' ■ (2.28) 

Here and in the following we append to the quasinormal frequencies the angular indices (I, m) and a third index 
(the "overtone" index N = 0, 1, 2 . . . ) that sorts frequencies by the magnitude of their imaginary part: in the limit 
a/M = 0, modes with larger iV have larger imaginary part and damp faster. The angular momentum-dependent 
quantity a m WNN'( a ) will be determined below, and the Kronecker symbol 5 m ,m' comes from the e lm( ^-dependence of 
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FIG. 7: Modulus | s £fj m | and phase s <t>im of the angular eigenfunction for the fundamental mode (thick lines) and first overtone 
(thin lines). The black hole's angular momentum a/M = (black, solid lines), 0.8 (red, dashed lines) and 0.98 (green, dot- 
dashed lines). The phase (hence the imaginary part) usually grows with the imaginary part of the mode frequencies, so it is 
larger for higher overtones. The modulus of the eigenfunction gives an intuitive picture of the angular pattern of the radiation: 
the radiation is almost entirely concentrated around x — cos# = 1 (-1) for corotating (counterrotating) modes, respectively. 
For axisymmetric perturbations (m = 0) most of the radiation is on the equatorial plane, x = cos 9 = 0. 
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the harmonics. For a/M — the spin- weighted spheroidal harmonics reduce to spin- weighted spherical harmonics 

s Y lm 0, for which 

-2^m -2Yi' m 'dQ. = 5i ; i<5 m , m i . (2.29) 

We evaluated SWSHs and their scalar products by two independent, approximate calculations. 

In the first calculation we use the sixth-order polynomial approximation for the separation constant s A\ m computed 
by Press and Teukolsky (Table I in Ref. 36]). We plug this expansion in the spheroidal wave equation and integrate 
the equation numerically imposing the appropriate boundary conditions at one extremum of the integration interval. 
In principle we should expect the agreement with Leaver's method to be good only for small values of au. Surprisingly, 
results for the eigenfunctions and the scalar products turn out to be very close also for large aw, where the polynomial 
approximation of the "true" separation constant is not so accurate. Even for a/M — 0.98 the approximate scalar 
products agree with Leaver's method within ~ 5%. 

As a second independent, analytic check we solved the spheroidal wave equation expanding in powers of au> and 
using standard perturbation theory. For auj = the solutions are the spin- weighted spherical harmonics |12|. The 
next order correction was found by Press and Teukolsky (Eq. (3.7) of Ref. [3(|). Using their approximate solution we 
can show that, to leading order, 

a m wNN'(a) ^ 4a [u* mN 0{l,l' ,m) + ui> m ' N >9(l' ,l,m)] , , (2.30) 

where ot m ii'NN'(fl) has been defined in Eq. (|2.28|l . 



9(1,1', m) = V2 jJ_l (/lmOU'm)(a20U'2) , (2.31) 

and (jij2?7T.im2 1 JM) is a Clebsch-Gordan coefficient. For I — V the same first-order perturbative analysis gives 

Oimll'nn'W — 1. 

In Table 0] we compare values of the scalar products obtained using our two approximate methods with the results 
of Leaver's method. There we focus on the worst possible scenario of a nearly maximally rotating Kerr black hole 
(a/M = 0.98). Both approximations work really well, providing a powerful consistency check of the numerical 
calculation. 



TABLE I: Comparison between different methods to evaluate the scalar product of spheroidal harmonics for a nearly-maximally 
rotating Kerr black hole (a/M — 0.98). We computed all numbers for N = N' = 0, so we only list the corresponding values 
of I, I' and m. The most accurate result comes from Leaver's method, but an expansion of the separation constant in powers 
of aui, followed by a direct integration of the angular equation, is accurate to within a few percent. The first order analytic 
expansion lj2~3?7)l also gives a surprisingly good approximation to the numerical results. 



Scalar product 


(1,1', m) 


Leaver 


Polynomial 


Analytic 


(2,3, 2) 
(2,3,1) 
(2,3,0) 
(2,3,-1) 
(2,3,-2) 


0.0275 - 0.0144i 
0.0525 - 0.0258i 
0.0608 - 0.036L 
0.0547 - 0.0378i 
0.0398 - 0.0308i 


0.0272 - 0.0150i 
0.0517 - 0.0273i 
0.0598 - 0.0379i 
0.0539 - 0.039L 
0.0394 - 0.0314i 


0.0279 - 0.0154i 
0.0510 - 0.0269i 
0.0591 - 0.0375i 
0.0543 - 0.0394i 
0.0409 - 0.0327i 


(3,4, 2) 
(3,4,1) 
(3,4,0) 
(3,4,-1) 
(3,4,-2) 


0.0393 - 0.0204i 
0.0487 - 0.0289i 
0.0487 - 0.0335i 
0.0440 - 0.0338i 
0.0363 - 0.0303i 


0.0389 - 0.0212i 
0.0482 - 0.0298i 
0.0484 - 0.0342i 
0.0438 - 0.034L 
0.0363 - 0.0304i 


0.0345 - 0.0187i 
0.0452 - 0.0279i 
0.0480 - 0.0338i 
0.0457 - 0.0355i 
0.0395 - 0.033L 



For reference, the constants otmll'nn' (a) computed using Leaver's method are given in Tables [n] and Ifffl for fixed 
values of m = m' (for m ^ ml the scalar product is zero). In the context of gravitational wave data analysis, it is 
most useful to consider scalar products between the quasinormal modes which are more likely to be excited during 
collapse or merger: say, the first three Kerr overtones (N, N') = 0, 1, 2 with (I, I') = 2, 3, 4 and \m\ < 2. Tables ITT1 
and lllll show that, to a reasonable level of approximation, we can assume a m u> nn r(a) ~ 8u> , that is, 

/ -2>S)m* {aUlmN)-2Sl'm' (O^l'm'N' )<m ~ Si (2.32) 

The errors on this approximate relation increase with the black hole's angular momentum a/M. They are typically 
larger for small m, but in most cases the approximate formula is valid to within ~ 10%, at least for cases of physical 
interest (say, for I < 5). 
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TABLE II: Scalar products of the spin weighted spheroidal harmonics for a/M — 0.80. To save space, we omit the leading 
zeros and the i's in the imaginary parts (so, for example, .9994 ± .0255 actually means 0.9994 ± 0.0255i). Entries replaced by 
an asterisk can be obtained from the symmetric entries by complex conjugation: i.e., a m ui NN i(a) = ct* ml , lN , N (a). 



m 


a/M = 0.80 


2 


Z' = 2 


f = 3 


I' =4 


I 


JV 


JV = 





N' = 1 


JV' = 2 


JV' = 


JV' = 1 


N' = 2 


JV' = 


JV' = 1 


JV' = 2 







1 




.9994-.0255 


.9976-.0510 


.0339-. 0221 


.0359-.0443 


.0392-.0665 


-.0016-.0005 


-.0024-.0007 


-.0036-.0011 


2 


1 






1 


.9994-.0255 


.0370-. 0441 


.0405-.0661 


.0454-.0880 


-.0026-.0013 


-.0039-.0018 


-.0057-.0024 




2 






* 


1 


.0422-. 0665 


.0473-.0881 


.0538-. 1096 


-.0044-. 0024 


-.0062-.0032 


-.0085-.0043 











* 


* 


1 


.9984+.0432 


.9937+.0869 


.0357-.0260 


.0389-.0503 


.0449-. 0743 


3 


1 


* 




* 


* 




1 


.9984+.0439 


.0352-.0535 


.0394-. 0778 


.0465-. 1016 




2 






* 


* 






1 


.0335- .0811 


.0388-. 1054 


.0469-. 1291 











* 


* 




* 


* 


1 


.9960+.0809 


.9842+. 1614 


4 


1 






* 


* 


* 


* 


* 


* 


1 


.9960+.0813 




2 






* 


* 




* 


* 


* 


* 


1 


1 


V = 2 


V = 3 


I' = 4 


/ 


JV 


N' = 





JV' = 1 


JV' = 2 


N' = 


N' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 







1 




.9995+.0074 


.9980+.0158 


.0455-.0295 


.0485-.0581 


.0540-.0868 


-.0023-.0007 


-.0033-.0010 


-.0049-.0015 


2 


1 






1 


.9995+.0083 


.0484-.0600 


.0533-.0883 


.0607-.1165 


-.0038-.0020 


-.0056-.0026 


-.0079-.0035 




2 






* 


1 


.0528-.0913 


.0596-.1193 


.0688-. 1468 


-.0066-.0038 


-.0091-.0048 


-.0121-. 0062 







* 




* 


* 


1 


.9968+.0677 .9869+. 1356 


.0390-.0291 


.0432-.0554 


.0512-. 0811 


3 


1 


* 




* 


* 


* 


1 


.9967+.0686 


.0372-.0608 


.0424-. 0871 


.0514-. 1127 




2 


* 




* 


* 


* 




1 


.0328-.0927 


.0390-. 1191 


.0491-. 1448 







* 




* 


* 






* 


1 


.9944+.0978 


.9776+. 1949 


4 


1 


* 




* 


* 


* 


* 


* 


* 


1 


.9943+.0985 




2 


* 




* 


* 




* 


* 


* 


* 


1 





I' = 2 


f = 3 


I' =4 


I 


JV 


N = 





N' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 







1 




.9985+.0420 


.9937+.0870 


.0484-. 0329 


.0524-.0635 


.0603-.0940 


-.0025-.0007 


-.0035-.0009 


-.0049-.0013 


2 


1 






1 


.9983+.0451 


.0504-.0690 


.0563-.0992 


.0661-. 1291 


-.0045-. 0022 


-.0062-.0028 


-.0085-.0037 




2 






* 


1 


.0520-. 1074 


.0599-. 1372 


.0717-. 1665 


-.0082-. 0044 


-.0109-.0055 


-.0139-.0070 











* 


* 


1 


.9946+.0924 


.9782+. 1853 


.0386-.0300 


.0436-.0562 


.0535-.0817 


3 


1 


* 




* 


* 




1 


.9944+.0942 


.0353-.0639 


.0413-.0902 


.0522-. 1158 




2 






* 


* 






1 


.0276-.0981 


.0346-. 1247 


.0465-. 1506 











* 


* 




* 


* 


1 


.9927+. 1141 


.9704+.2275 


4 


1 


* 




* 


* 




* 


* 


* 


1 


.9925+.1154 




2 


* 




* 


* 




* 


* 


* 


* 


1 




-1 


r = 2 


/' =3 


I' = 4 


I 


JV 


JV = 





N' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 







1 




.9963+.0789 


.9834+. 1658 


.0439-.0321 


.0488-.0606 


.0590-.0892 


-.0023-.0005 


-.0030-.0005 


-.0040-.0007 


2 


1 






1 


.9954+.0877 


.0444-. 0695 


.0511-.0977 


.0630-. 1257 


-.0043-.0020 


-.0058-.0023 


-.0075-.0031 




2 






* 


1 


.0410-. 1117 


.0497-. 1397 


.0637-. 1673 


-.0087-.0040 


-.0110-.0050 


-.0134-.0065 







* 




* 


* 


1 


.9922+. 1166 


.9676+.2350 


.0354-.0289 


.0410-.0533 


.0522-.0769 


3 


1 


* 




* 


* 


* 


1 


.9916+. 1208 


.0308-.0625 


.0373-.0872 


.0495-. 1112 




2 


* 




* 


* 


* 




1 


.0199-.0966 


.0274-. 1220 


.0406-. 1468 







* 




* 


* 




* 


* 


1 


.9908+.1294 


.9628+.2586 


4 


1 


* 




* 


* 




* 


* 


* 


1 


.9905+.1320 




2 


* 




* 


* 






* 


* 


* 


1 




-2 


I' = 2 


I' = 3 


l'=4 


I 


JV 


N' = 





JV' = 1 


JV' = 2 


JV' = 


N' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 







1 




.9927+. 1170 


.9656+.2523 


.0327-.0257 


.0376-.0476 


.0486-.0700 


-.0017-.0002 


-.0020-.0001 


-.0025-.0002 


2 


1 






1 


.9899+.1377 


.0319-.0577 


.0382-.0794 


.0505-.1016 


-.0032-.0013 


-.0041-.0015 


-.0050-.0020 




2 


* 




* 


1 


.0237-.0963 


.0319-.1183 


.0461-. 1408 


-.0071-.0027 


-.0086-.0034 


-.0101-. 0046 







* 




* 


* 


1 


.9895+.1397 


.9556+.2841 


.0298-.0255 


.0354-.0465 


.0469-.0668 


3 


1 


* 




* 


* 




1 


.9882+. 1479 


.0245-.0560 


.0310-.0775 


.0434-. 0984 




2 


* 




* 


* 




* 


1 


.0118- .0872 


.0190-. 1097 


.0324-. 1318 







* 




* 


* 




* 


* 


1 


.9890+.1436 


.9551+.2879 


4 


1 


* 




* 


* 




* 


* 


* 


1 


.9884+. 1478 




2 


* 




* 


* 




* 


* 


* 


* 


1 
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TABLE III: Scalar products of the spin weighted spheroidal harmonics for a/M — .98. To save space, we omit the leading 
zeros and the i's in the imaginary parts (so, for example, .9997 ± .0213 actually means 0.9997 ± 0.0213i). Entries replaced by 
an asterisk can be obtained from the symmetric entries by complex conjugation: i.e., a m ui NN i(a) = ct* ml , lN , N (a). 



m 












a/M = 


.98 








2 


l' = 2 


J = 3 


J' =4 


I 


jV 


JV = 





JV' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 







1 




.9997-.0213 


.9987-.0427 


.0275-.0144 


.0277-. 0294 


.0280-.0437 


-.0027+.0000 


-.0033+.0001 


-.0041+.0004 


2 


1 


* 




1 


.9997-. 0214 


.0284-.0277 


.0293-.0428 


.0303-.0570 


-.0032-.0001 


-.0040-.0000 


-.0052+.0001 




2 








1 


.0299-.0411 


.0315-.0560 


.0331-.0702 


-.0039-.0003 


-.0051-.0003 


-.0066-.0002 













* 


1 


.9993+.0269 


.9975+.0524 


.0393-.0204 


.0399-.0412 


.0402-.0613 


3 


1 






* 


* 


* 


1 


.9994+.0256 


.0400-.0399 


.0413-.0605 


.0422-.0804 




2 








* 




* 


1 


.0407-.0583 


.0425-.0787 


.0440-.0984 













* 




* 


* 


1 


.9973+.0662 


.9897+. 1292 


4 


1 






* 


* 


* 


* 


* 


* 


1 


.9975+.0633 




2 






* 


* 


* 


* 


* 


* 


* 


1 


1 


I' = 2 


Z' = 3 


V = 4 


I 


jV 


JV' = 





JV' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 







1 




.9997+.0031 


.9990+.0063 


.0525-.0258 


.0532-.0528 


.0533-.0802 


-.0036-.0006 


-.0047-. 0009 


-.0064-.0013 


2 


1 






1 


.9998+.0031 


.0569-.0492 


.0590-.0758 


.0605-. 1028 


-.0046-.0017 


-.0062-.0023 


-.0085-.0029 




2 






* 


1 


.0638-.0685 


.0669-.0947 


.0695-.1210 


-.0057-.0031 


-.0078-. 0040 


-.0105-.0050 











* 


* 


1 


.9971+.0632 


.9886+. 1255 


.0487-.0289 


.0518-.0553 


.0574-.0814 


3 


1 






* 


* 


* 


1 


.9971+.0626 


.0489-.0601 


.0531-.0862 


.0598-. 1119 




2 








* 




* 


1 


.0492-.0915 


.0546-. 1172 


.0623-. 1424 













* 


* 


* 


* 


1 


.9942+.0996 


.9768+. 1976 


4 


1 






* 


* 




* 


* 


* 


1 


.9942+.0994 




2 


* 




* 


* 




* 


* 


* 


* 


1 





I' = 2 


Z = 3 


Z = 4 


/ 


jV 


AT = 





N' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 







1 




.9983+.0446 


.9927+.0918 


.0608-.0361 


.0653-.0690 


.0734-. 1017 


-.0039-.0010 


-.0050-.0014 


-.0069-.0020 


2 


1 


* 




1 


.9981+.0474 


.0639-.0758 


.0708-. 1081 


.0812-.1399 


-.0061-.0031 


-.0083-.0039 


-.0110-.0052 




2 








1 


.0683-. 1179 


.0776-. 1493 


.0905-. 1798 


-.0103-.0063 


-.0134-.0079 


-.0169-.0100 











* 


* 


1 


.9934+.1025 


.9733+.2043 


.0487-.0335 


.0547-.0616 


.0660-.0886 


3 


1 






* 


* 


* 


1 


.9932+. 1035 


.0449-.0723 


.0522-. 1005 


.0649-. 1275 




2 


* 






* 




* 


1 


.0361-. 1111 


.0449-. 1396 


.0590-. 1670 











* 


* 




* 


* 


1 


.9907+.1282 


.9628+.2543 


4 


1 








* 


* 


* 


* 


* 


1 


.9906+.1289 




2 






* 


* 


* 


* 


* 


* 


* 


1 




-1 


/' = 2 


r = 3 


V = 4 


/ 


TV 


JV' = 





JV = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 







1 




.9948+.0935 


.9769+. 1954 


.0547-.0378 


.0618-.0697 


.0760-. 1007 


-.0036-.0007 


-.0045-.0008 


-.0058-.0012 


2 


1 


* 




1 


.9936+. 1032 


.0544-.0832 


.0641-. 1146 


.0808-. 1450 


-.0064-.0028 


-.0083-.0035 


-.0104-. 0047 




2 






* 


1 


.0481-. 1339 


.0609-. 1652 


.0807-. 1951 


-.0126-.0056 


-.0155-.0072 


-.0185-.0096 











* 


* 


1 


.9891+.1376 


.9553+. 2752 


.0440-.0338 


.0517-.0606 


.0667-.0857 


3 


1 






* 


* 




1 


.9884+. 1414 


.0370-. 0742 


.0462-. 1016 


.0628-. 1275 




2 








* 




* 


1 


.0212-. 1144 


.0318-. 1431 


.0501-. 1703 













* 




* 


* 


1 


.9874+.1519 


.9491+.3017 


4 


1 


* 




* 


* 


* 


* 


* 


* 


1 


.9870+. 1542 




2 








* 




* 


* 


* 


* 


1 




-2 


I' = 2 


Z = 3 


Z = 4 


/ 


JV 


JV' = 





JV' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 


JV' = 


JV' = 1 


JV' = 2 







1 




.9891+.1430 


.9478+.3089 


.0398-.0308 


.0475-.0554 


.0641-.0797 


-.0025-.0003 


-.0030-.0001 


-.0035-.0003 


2 


1 


* 




1 


.9844+. 1704 


.0374-.0709 


.0473-. 0955 


.0662-. 1196 


-.0049-.0018 


-.0060-. 0021 


-.0071-.0031 




2 






* 


1 


.0223-. 1195 


.0353-. 1452 


.0575-. 1705 


-.0110-.0036 


-.0129-. 0048 


-.0147-.0070 











* 


* 


1 


.9847+.1686 


.9355+.3406 


.0363-.0303 


.0445-.0535 


.0609-.0748 


3 


1 


* 






* 




1 


.9828+. 1783 


.0280-.0676 


.0376-.0918 


.0556-. 1145 




2 






* 


* 




* 


1 


.0082-. 1048 


.0191-. 1310 


.0388-. 1561 











* 


* 




* 


* 


1 


.9843+.1718 


.9359+.3423 


4 


1 






* 


* 


* 


* 


* 


* 


1 


.9834+.1766 




2 






* 


* 




* 


* 


* 


* 


1 
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III. HIGHER-DIMENSIONAL SPHEROIDAL HARMONICS 



A. Series solution 



Scalar (n + 4)-dimensional spheroidal harmonics are defined as solutions of the equation 



sin 6 cos" 9 \d9 



,dS; 



kjm 



de 



c cos 



'■8 -]{]+■ 



1) sec 2 9 + A 



kjm 



S k 



jm 



0. 



(3.1) 



In four dimensions there is only one possible rotation axis for a cylindrically symmetric space-time, so there is 
only one angular momentum parameter. In higher dimensions there are several possible choices of rotation axes. 
Therefore higher-dimensional black holes will have, in general, several angular momentum parameters, each referring 
to a particular rotation axis p9|. Here we focus on the simple case where there is only one rotation axis (and 
correspondingly, only one angular momentum direction). The angular separation constant A k j m now depends on 
three indices. The integer m = 0, ±1, ± 2 . . . comes from separation of the angle describing the azimuthal 
dependence of the perturbations around the symmetry axis. Since (|3.1[) is invariant under the replacement m — > — m 
we will restrict our calculations to m > 0. The number of extra dimensions n and the parameter j = 0, 1, 2, . . . 
are related with the eigenvalues of hyperspherical harmonics on the n-sphere, which are given by + n — 1). 

The parameter k(= 0, 1, 2 . . . ) plays the same role as r/2 in (|2.11() in the four-dimensional case, i.e. it labels the 
eigenvalues for fixed values of j and m. Eq. 1)3.1(1 describes the angular dependence of a scalar field perturbing a 
(n + 4)-dimensional Myers- Perry black hole with a single rotation parameter, and is a simple generalization of the 
four-dimensional case f30L l3ll 1321 133| . When n = 1 it also describes scalar perturbations of smooth geometries in the 
D1-D5 stringy system |35j . 

In the region of physical interest 9 £ [0, ir] the differential equation (|3.1|) has three regular singular points: 9 = 0, ir/2 
and 7r. We make the following change in the angular function: 



">kjm 



(sin6>) feo (cos6>) fel Vkjr 



(3.2) 



and using a Frobenius expansion for ykjm around 9 = 0, ir/2 and ir we find the characteristic exponents: fco = \m\, 
k\ = j or ki = 1 — n — j. The solution with k± = 1 — n — j is irregular at x = in the higher-dimensional case, so we 
take k\ = j. Inserting the expansion l|3.2(l into the differential equation l|3.1|l with ko = \m\ and k\ = j, we find that 
the function ykjm satisfies the differential equation 



(1 



dx 2 



(n + 2j){l-x 2 ) 



2(|m| + l)x 



dx 



(3.3) 



where S kjm = A k]m ~ j(j + n + 1) - \m\(n + \m\ + 2j + 1). 
We assume a series expansion for ykjm such that 



Sk 3 m = {sm9)\ m \cos9) j ^h p (cos 2 9) 

p=0 



I' 



(3.4) 



This series (if convergent) automatically satisfies the regul arity boundary conditions at 9 = 0, tt/2, tt. Substituting 

term recursion reli 

a di + (3 a = . 



arity 

(|3.4|) into H3.1fl we obtain the three-term recursion relation 33] 

a di + (3 a = , 

a p d p+1 + flpCLp + %d p -i = , (p = 1, 2, • • • ) 



where 

& P = -2(p+l)(2j + n + 2p+l) , 

h = (j + \m\+2p){j + n+\m\+2p+l)-A kjm , 

% = -c 2 - 

The continued fraction equation for the separation constant has the same form as Eq. 1)2. 9p : 

a 7i ai72 0:273 



(3.5) 



(3.6) 



/So 



Pi- P2- P 3 - 



0. 



(3.7) 
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B. Small-c expansion 



In this Section we provide an expansion of the eigenvalue Akj m in powers of c around c = 0. For c = 0, the 
eigenvalue A^j m is explicitly determined from the requirement that the series expansion has a finite number of terms, 
since otherwise it is divergent [l^. Imposing fa = for some integer k > we have, for c = (an alternative 
derivation is presented in Appendix [Al . 



If we identify 



A kjm = (2k + j + \m\)(2k +j + |m| + n + 1) . 



2k = l-(j+ \m\) . 



(3.8) 



(3.9) 



we can write Akjm = 1(1 + n + 1). The integer / is constrained to satisfy the condition I > (j + \m\), which is a simple 
generalization of the four-dimensional case. An important difference from the four-dimensional case is that now I 
cannot be any positive integer: for k to be a positive or zero integer, only even (odd) values of I are admissible when 
(j + \m\) is even (odd). 

For finite but small c we proceed as in Sec. Ill CI When c = the recursion relation l|3.7l) has a finite number of 
terms k. Therefore, expanding around c = we find it convenient (although strictly not necessary) to use the fc-th 
inversion of Eq. I|3.7|l : 



fa 



Oik-ilk ot k -2lk-i «o7i afc7fc+i afc+i7fc+2 



fa-x Pk-2- fa 
Now we expand the separation constant as a Taylor series: 



fa+l~ Pk+2~ 



(3.10) 



where 



.4 



kjm 



p=0 



fa 

h 
h 
h 
h 



l(l + n+l), 
0, 

-1 + 2j(j - 1) + 21(1 + 1) - 2m 2 + 2n(j + I) + n 2 



(21 + n- l)(2l + 3 + n) 



0, 

(j-l+\m\)(j + l-\m\+n-l) 



-1 + 2j(j - 1) + 21(1 + 1) - 2m 2 + 2n(j + I) + n 2 



16(21 + n - 3) (21 

4(21 + n- 3) 
y ' (2l + n-l)(2l + n + 3) 

(j - I + \m\ - 2)(j + l + n-\m\ + l) 



(2 + j-l + \m\)(l+j - \m\ + n-3) - 



16(2Z 



n)(2l 



■3) 2 



(j-l + \m\ - 4)(Z + j + n - |m| + 3) + 



h 



4(21 
0. 



-1 + 2j(j - 1) + 21(1 + 1) - 2m 2 + 2n(l + j) + n 2 
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l)(2l 



3) 



(3.11) 



(3.12a) 
(3.12b) 

(3.12c) 

(3.12d) 



(3.12e) 
(3.12f) 



The coefficients /i , . . . , /s were obtained substituting (|3.11() into Eq. I|3.10|l and expanding the resulting expression in 
powers of c. Higher order coefficients can be obtained easily, but the expressions are too lengthy to reproduce them 
here. The four-dimensional results can be obtained by setting j = n = 0. They are in agreement with j3]],[3|l and 
with our Eqs. I|2.16ajl - P35g| |. 



C. Large and real c (oblate case) 

We now proceed to calculate the asymptotic behaviour of the oblate HSHs (i.e., Eq. (|3.1|) with c € K) for large c, 
following the same method as in Sec. Ill Dl In order to find asymptotic solutions valid near the end-points (i.e. "inner 
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solutions") we make the change of independent variable u = 2c(l — x). The ordinary differential equation (|3.3(l then 
becomes 



c 



1 / 2 d2 Vkjm ,„,, | . n dy k jm u 2 \ (n + 2j)u f U\dy k j m 

Tc \ u ~W + 2(|m| + 1)u — ~ T Vk ^j - l2^u) V - Tc) — = ° 



(3.13) 



This equation is the same as the corresponding equation in four dimensions where 8 k j m plays the role of s Ai m , except 
for the presence of the last term in Eq. i|3.13[) . This last term can be neglected to leading order in c. The asymptotic 
solution of Eq. (|3.13|l valid near the end-points is therefore the same as the inner solution in four dimensions (given 
by Eq. (3.10) of 27] with h = 0), so that S kjm = -c 2 + 2qc + 0(1) = A kjm + 0(1) [cf. Eq. f2~T7jl ], 

Our task is now to determine the parameter q. We can find an asymptotic solution which is valid far from the 
end-points (the "outer" solution) as in [g^. The result is: 

S° k ? m (.x) = (1 - xT 1 / 2 ^"/ 2 [a kjm (l - xf' 2 (l + x)-^ 2 e +cx + b kjm (l - x)^l 2 {l + x)+^ 2 e~ cx ] (3.14) 

where a k jm and b k ±m are constants of integration. This solution is similar to the outer solution in four dimensions, 
i.e., Eq. (3.26) of 27] with h = 0. Unlike the four-dimensional case, l|3.14|l is not valid very near the origin x = 0. 
Despite this inconvenience, the matching of the inner and outer solutions in the overlap region proceeds exactly as in 
four dimensions: the number of extra dimensions n plays no part in the matching. 

The next step is to find the number of zeros of the HSHs. Equating this number with the number of zeros of the 
asymptotic solution determines the value of q. Refs. 0, @ show that scalar spheroidal harmonics (a special limit of 
HSHs with s = n = j = 0) are entire functions of c, so their number of zeros does not vary with c. In particular, 
their number of zeros is the same as that of the spherical harmonics (s = n = j = c = 0), which is known. Similarly, 
Refs. I45L 0] show that SWSHs (n = j = 0) are entire functions of c. Therefore their number of zeros is easily 
obtained, since it does not vary with c. The argument used in both cases relies crucially on two properties: 

(PI) the fact that the indicial equation at any possible regular singular point of the differential equation does not 
involve c, 

(P2) the continuity in x of the coefficients in the differential equation when the coefficient of the second order derivative 
is equal to one. 

We must be careful if we wish to use a similar argument to show that the number of zeros of the HSHs does not 
vary with c. First of all, the characteristic exponents <So and ct\ for the HSH equation do not depend on c, and the 
differential equation (|3.3I) does not possess any singular points. Therefore the differential equations for S k j m an d ykjm 
both satisfy the property (PI) above. However, the coefficients of either differential equation are not continuous at 
x = when the coefficient of the second order derivative is equal to one. Therefore, property (P2) is not satisfied 
by either differential equation. Nevertheless, we can still apply the argument separately to the regions x E (—1,0) 
and x G (0, +1), where both properties (and other necessary ones) are satisfied by both differential equations. Note 
that no zero of ykjm can "move" from either region x G (—1,0) or x G (0, +1) to the other as c varies, because this 
function cannot possess a zero at x — 0. We then conclude that the number of zeros of ykjm in either interval (—1, 0) 
or (0, +1) [and therefore in the whole region (—1, +1)] is independent of c. In particular, it is equal to the number of 
zeros of ykjm in the corresponding region for c = 0. So the number of zeros of the HSHs for x G (—1, +1) for any real 
c can be obtained from the special case c = 0. 

In Appendix [X] we find an exact solution of the HSH equation for c = and show that the number of zeros of the 
HSHs in x G (—1, +1), excluding x — 0, is equal to 2k. The HSHs possess a zero at x — if j > and do not if j = 0; 
however, we will not need this property. As a check of the analytic results we performed numerical calculations of the 
eigenfunctions using Eq. I|3.4f> : some of these numerically computed, normalized eigcnfunctions are shown in Fig. EI 

The matching of the number of zeros of the HSH with that of the asymptotic solution proceeds now as in four 
dimensions. Because of the forms of the inner and outer asymptotic solutions, the number of zeros of the asymptotic 
solution (excluding x = 0) has the dependence on q given in Eq. (4.4) of |27J with h = (h in that paper is what 
we call s here) . We can equate this number to the number of zeros of the HSH (excluding possible zeros at x = 0), 
which is 2k = I — \m\ — j. 

The resulting equation is the same as Eq. (4.4) in |2JJ with h = 0, and with zq replaced by j. The value of the 
parameter q is therefore 



q = l-j + l = 2k+ \m\ + 1 



(3.15) 
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FIG. 8: Some representative oblate eigenfunctions Skjm for different values of c. The plot on the left refers to n = 2. For c — 
the eigenfunction can be obtained analytically (see Appendix lAl with the result Si 21 = [45(1 — x 2 )/32] (7 — tlx 2 ) x 2 . In 
the plot on the right n = 1, and for c = we have S201 = [693(1 - x 2 )/380] 1/2 (l - 8a; 2 + lOx 4 ) . 

The asymptotic behaviour of the separation variable Akj m is then given by 

A kjm = -c 2 + 2(1 - 3 + l)c + 0(1) . (3.16) 

There is an imprint of the extra dimensions (through the parameter j) already at order O(c). Notice that the above 
asymptotic expansion (|3.16|) does not agree in the limit n — j — with the four-dimensional result, Eqs. Q2.17jl and 
(|2.19() with s — 0. The reason is that the present analysis is based on picking the characteristic exponent k\ — j while 
disregarding fcx = 1 — n — j in Eq. I|3.2|) . The latter value should not be disregarded when n = j = 0, in which case 
ai can take on both values, and 1. 




FIG. 9: Oblate eigenvalues of the HSH with k = n = 1 for m = 1 (left), m = 2 (right), and different values of j. 

To conclude this Section, we back up our analytic results by numerical calculations of the oblate eigenvalues 
obtained solving Eq. i|3.7fl . In Fig. |5]we plot the eigenvalues for n = k = 1, m = 1, 2 and j = 0, 1, 2; these values 
of (n, k,j, m) have been chosen to match those in Figs. 2-7 of [3l|. We checked that the asymptotic behavior of the 
oblate eigenvalues, both in the limit c — > and in the limit c — > 00, is in agreement with the analytic predictions. 
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D. Large and pure-imaginary c (prolate case) 



In this Section we complete the picture computing a series expansions of the eigenvalue Akj m for large and pure- 
imaginary values of c. The present analysis is similar to the one in Sec. Ill El 

We first find an approximation valid near the origin when |c/| — > oo by introducing a new independent variable 
u = y/2 1 cj | x and using again the transformation (|3.2|1 (with fco = \m\ and k\ = j). The resulting ordinary differential 
equation for the function ykjm is the same as Eq. Ij3.3l) . but in terms of the new variable u: 



2 



2\c I \{n + 2j) 



- (n + 2j + 2\m\ + 2) u 



2 



d i X \ Cl \ U K n r> 

j- + Ofejm ^ — f = ' \ / 



Now we apply the same approximation as in Sec. Ill El In the limit |c/| — * oo and in the inner region l/|c/| -C |x| <C 1, 
Eq. I|3.17[l reduces to 



d 2 (n + 2j) d u 2 



fe = 0, (3.18) 



where we have set Skjm ~ 2g|cj| + (|c/|°), and g is an unknown parameter that we shall now determine, not to be 
confused with the parameter introduced in (|3.15|) . With respect to the four-dimensional case, Eq. I|3.18|) contains an 
extra term with a first derivative. The presence of this term means that the inner solution is not a parabolic 
cylinder function, as in the four-dimensional case. By applying the change of variables 

y^ m {u)=u-^/ 2 v{u), (3.19) 

Eq. p.lSjl becomes 

id 2 (n + 2j)(n + 2j-2) u 2 \ 

\d^ + q An 2 Tf = °- (3 ' 20) 

The solution of this equation can be expressed in terms of a Laguerre polynomial L*~** ^tJ. It can be checked that 
the asymptotic approximation of the prolate HSHs valid in the inner region l/|c/| <C |af| <C 1 as |cj| — > oo, which is 
regular at x — 0, is 



^(u)= e -» S / 4 4 5) (u 2 /2) , (3.21) 



where 



n + 2J-1 2g - (n + 2j + 1) 
a = and n = . (3.22) 

The number of positive zeros of the Laguerre polynomial Lj*\x) is equal to n if (as in our case) a > — 1 0. ^ can 
easily be checked that all these zeros lie within the inner region. The number of zeros of y™^ m in the inner region [or 
equivalently in the region x € (— 1, +1)] is therefore equal to 2n. 

For \a\ — ► oo, an outer solution which is a valid approximation of the prolate HSH near the end-points can be 
obtained using a WKB-type expansion, as in and in Sec. Ill El The result is 



b k 



kjm 



(3 23) 

(i - x 2 y lfi x- 1 / 2 -^ (i + yr^jV^i^ 3 



where akjm and 6fcj m are constants of integration. For n = j = 0, (|3.23|) reduces to the corresponding outer solution 
(|2.2(il) in the prolate SWSH case if s = and we replace L + 1/2 there by q. Furthermore l|3.23|l is valid in the same 
region as l|2.26|l : for |x| >> 1/a/RzT and \x\ ~ 1 - <3(|c/| e ), with -1 < e < 0. Matching the outer solution Ij3.23jl with 
the inner solution in the overlap region we find 

J (a + n — z)!n!(a + n)\ V / 
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This outer solution has no zeros. Like in the prolate SWSH case with s = 0, all zeros of the asymptotic solution for 
\cj\ — ► oo are provided by the inner solution. We want to equate their number to the number of zeros of the exact 
solution ykjm in the region x £ (—1, +1), which is given in Appendix 1X1 for c = 0, is independent of c 2 £ (— oo, +oo) 
and is equal to 2k. Therefore 2q = 4fc + n + 2j + 1, and the separation constant for large |c/| is 

A kjm = 6 kjm + O (| CJ |°) - (4fc + n + 2j + l)\cj\ + O (| Cj |°) = [2(1 - |m|) + n + 1] \a\ + O (|c z [°) . (3.25) 




FIG. 10: Prolate eigenvalues of the HSH with k = n = 1 for m = 1 (left), m = 2 (right), and different values of j. 

The expansion (|3.25|l is in agreement with the eigenvalues obtained from a numerical solution of Eq. I|3.7() in the 
prolate case, some of which are plotted in Fig. 



IV. CONCLUSIONS 

In this paper we provided a complete numerical and analytic analysis of four-dimensional spin- weighted spheroidal 
harmonics (SWSHs) and of higher-dimensional spheroidal harmonics (HSHs). For reference we list below what we 
regard as the most useful results, be they known in the existing literature or originally derived in this paper. 

For SWSHs, the main results are: 

1) Eq. (|2.5() gives a series representation of the SWSHs, and a simple computational algorithm to determine 
eigenvalues and eigenfunctions by recursion; 

2) Eq. I|2.13|) is the series expansion of the eigenvalues for small c; 

3) Eq. I|2.17|l is the series expansion of the eigenvalues for large cgl (oblate case); 

4) Eq. H2.21fl (for spin 0) and Eq. (|2.24|) (for general s, but to lower order) are series expansions of the eigenvalues 
for large c when ic e M. (prolate case) ; 

5) Eq. H2.32fl and, more accurately, Eq. H2.30fl provide analytic approximations of the scalar products between 
SWSHs corresponding to slowly damped eigenfrequencies of Kerr black holes. 

For HSHs, the main results are: 

1) Eq. (|3.4|l gives a series representation of the HSHs, and a simple computational algorithm to determine eigen- 
values and eigenfunctions by recursion; 

2) Eq. I|3.11|) is the series expansion of the eigenvalues for small c; 
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3) Eq. 13.16|l is the series expansion of the eigenvalues for large c £ M. (oblate case); 

4) Eq. 13.25|) is the series expansion of the eigenvalues for large c when ic € R (prolate case); 

5) Eq. (|A12I) is an exact solution for the eigcnfunctions of the HSH equation for c = in terms of Jacobi polynomials; 
Appendix iBl contains another exact solution of the HSH equation for special values of the arguments. 

In our opinion, the main unsolved theoretical problem is to clarify why Eq. (|2.24|) is in good agreement with the 
numerics for \s\ ^ 0, even though the number of zeros of the eigenfunction seems to depend on c when c is not 
real. Given the excellent accuracy of Eq. I|2.30|l in the four-dimensional case, it would be desirable to extend the 
perturbative formalism developed by Press and Teukolsky [3G| to higher dimension. It could also be useful to push 
the expansions (|2.24|) . (|3.16|) and (|3.25|l to higher order, and to extend the present understanding of the location of 
the spheroidal branch points [l(| to the case s ^ 0. Finally, it would be an exciting and interesting challenge to 
extend our results for higher-dimensional scalar harmonics to higher-dimensional harmonics with s ^ 0. We leave 
these issues for future work. 
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APPENDIX A: EIGENVALUES, EIGENFUNCTIONS AND ZEROS OF THE HSH EQUATION FOR c = 

In the c = limit the HSH equation can be solved exactly, and the eigenvalues can be computed analytically. The 
general solution for the eigenvalues was first found by Ida et al. |3lj . However, their method relies on the continued 
fraction representation and does not provide an explicit analytic expression for the eigenfunctions. The purpose of 
this Appendix is to compute explicitly the eigenfunctions and eigenvalues for c = 0. 

Let us introduce a new wavefunction "f: 

S fcj>re = (sin0)H(cos0)W*. (Al) 
Setting c = and z = (cos#) 2 in l|3.1|l we get 

z(l - zR 2 * + ij + - + j + \m\j zj d z V - ((j + m)(j + |m| + n + l)/4 - A kjm /A) * = . (A2) 

This is an hypergeometric equation: 

z(l - z)dl^> + (7 - (a + (3 + l)z) d z ^> - (a/3) * = , (A3) 

with 



n + 1 ll. , , n+1 I , fn + 1 N 2 




7 = •?+— ' a = 2 y + |m| + ~~ ~ + y kjm 

The general solution for z e (0, 1) can be written as 

* = Az l ~"<F[a - 7 + 1, (3 - 7 + 1, 2 - 7, z] + BF[a, 0, 7, z] . (A5) 
For small z (sin 8 ~ 1) F ~ 1, and we get 

~ A cos 0- 2 ^ + B , S k jm ^ Acos8-^ + B . (A6) 
Regularity requires A = 0. To analyse the point z = 1 we use the property 
T7ir t 1 n Nc-a-b£[c]r^a + fe-c] T[c]T[c-a-b] 

F[a,b,c,z] = (l-z) rTTFTIi F[c-a,c-b,c-a-b+l,l-z]+—. t— jtF [a, b, a+b-c+1, 1-z] . A7 

1 [a\l [b\ I [c — a\L [c — b\ 
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So for z ~ 1 



m -2| H W« + M , r[ 7 ]r[ 7 - a - /3] 



r[«]r[/3] r[ 7 - a]r[ 7 -/3] 

and requiring regularity of the wavefunction we get a = — k or ft = — k (fe = 0,l,2,...). This leads to 

A kjm = (2k + j + \m\) (2k + j + |m| + n + 1) , (A9) 

in agreement with the result found by Frolov and Stojkovic for n — 1 30] and by Ida et al. for general n 31] . The 
eigenfunctions can be written in terms of hypergeometric functions: 

S kjm = (sm9)^(cos9)^F[a, /?, j,cos 2 9} , (A10) 

For reference we list some eigenfunctions for n = 1, obtained from 1A10|) and then normalized: 

S 00 o = (1/2) 1 / 2 , Soio = (3/2) ^ 2 cos 9, S QO i = (3/4) 1 / 2 sin9 , 

5on = (15/4) 1/2 cos(9sin6l, 5 02 i = (35/4) 1/2 cos 2 sin , 

Sioo = (15/14) 1/2 (1 - 2cos 2 0) , Swi = (21/16) 1/2 (l - 3 cos 2 6») sin 6> , 

S ln = (63/4) 1/2 cos(9(l-2cos 2 (9)sin6l, S ln = (31185/464) 1/2 cos 2 9 (3 - 5 cos 2 9) sin 6 , 

5io 2 = (315/208) 1/2 sin 2 6»(l-4cos 2 6i) , 5u 2 = (3465/592) 1/2 cos (2 - 5 cos 2 9) sin 2 9 , 

#122 = (45045/368) 1/2 cos 2 6»(l-2cos 2 6»)sin 2 6l, S 20 2 = (45 45/ 22 2 08) 1/2 [sin 2 9 (l - 10 cos 2 9 + 15 cos 4 9)] . 

To convert hypergeometric functions to Jacobi polynomials T"^^(z) we can use the property [43| : 

F[a, -n, c, z] = j^V^ 1 ■- c -»(l - 2z) , (All) 

where (c)„ = Y\c + n]/T[c] is the Pochhammer symbol. As a result, the (non normalized) wavefunctions can be 
expressed as 

S = (sin 9)^ (cos 9)^ V^- (1-2 cos 2 9), (A12) 

where a = (k — l)/2 + j and [3 = \m\. 

We can now use this analytic expression of the eigenfunctions to determine the number of zeros in the interval of 
physical interest. Define Klein's symbol 

( if u < 

E(u) = I [u] if u positive and nonintegral (A13) 
I u — 1 if u = 1,2,..., 

where [] denotes the floor function, and 

X(a ,0) = E Q(|2n + a + + 1\ - \a\ - |/3| + 1)) . (A14) 
Then the number of zeros N(a , (3) of V° '^(z) in the interval [—1,1] is given by 0] 

2[(X+l)/2] for (-l)' r 

N(a,f3) = { }J1 R { (A15) 

2[X/2]+l for (-1)' 

Specializing to our case, we get that the number of zeros of V£ ' (z) in the (—1,1) interval is k. But '^(1 — 2 cos 2 9) = 
Vj? '^(— cos 26*), so in the interval (0 , tt) we have 2k zeros. We know from Eq. I|3.2|l that yi m does not possess a zero 
at cos 9 — 0,±1. This is ratified here by the fact that (— cos 29) does not possess a zero at cos6> = nor at 
cos# = ±1 because neither a nor j3 can take on the values —1, —2, . . . , — k when k, j, \m\ > 0. 
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APPENDIX B: SOLUTION OF THE HSH EQUATION FOR SPECIAL VALUES OF THE ARGUMENTS 

It is well known that the four-dimensional scalar spheroidal harmonics admit an exact solution when m = 1 and 
oAu — — c 2 [lj. For these values the spheroidal equation reduces to 

<9 2 *-c 2 <& = G, (Bl) 

where \fr = (x 2 — l) 1 / 2 qSh- Then the solution is elementary: <!> = sin(icx), with 

ic = lTr, 1 = 0,1,2... (B2) 

In an attempt to generalize this simple particular solution to the (n + 4)-dimensional case, we introduce a wavefunction 



V = x n/ Wx 2 ~lS k3m . (B3) 

We find that ^ admits the same elementary solution as in the four-dimensional case only if j = 0, m = 1, n = 2, 
Akoi = — c 2 — 2, and then c takes the value (|B2J| . We could not find other values of the arguments leading to such a 
simple solution. 
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